Comment on "Macrospopic Equation for 
the Roughness of Growing Interfaces in 
Quenched Disorder" 

In a recent Letter jl]] Braunstein and Buceta intro- 
duced a 'macroscopic' equation for the time evolution of 
the width of interfaces belonging to the directed perco- 
lation depinning (DPD) universality class @. From nu- 
merical simulations of the DPD model, they inferred an 
ansatz (Eq.(l) in Ref. for the time derivative of the 
interface width (called DSIW in Ref. 0) at the depin- 
ning transition. Braunstein and Buceta found that their 
formula fitted the numerical data at the depinning tra- 
sition, for q c = 0.539 and (3 — 0.63, with the appropriate 
election of some arbitrary constants. 

Here we argue that, contrary to what it is claimed in 
Ref. 0, Braunstein and Buceta's formula does not de- 
scribe the 'macroscopic' behaviour of the interface. The 
formula proposed in Ref for the DSIW is an approxi- 
mation to the very short times regime (when less than one 
layer has been completed) , which is not significant for the 
description of the surface dynamics at large scales. We 
obtain analitically the short time behaviour of the DPD 
model, which is valid for any q and explains the apper- 
ance of an exponential term in the formula of Ref. for 
the DSIW. 

Let us consider the DPD model in a system of size 
L and a density q of blocked cells (p = 1 — q den- 
sity of free cells). We are interested in the very short 
times regime when the first monolayer still has not been 
completed, i.e. the number of growth attempts N is 
N <C L (this corresponds to times t = N/L « 1). In this 
regime, the probability of having a column i with height 
hi > min(hi-i, ^i+i)+2 is negligible and the columns are 
growing almost independently. The growth at this early 
stage can be seen as a random deposition (RD) process 
[S in which every column grows in one unit with prob- 
ability p/L. The short time regime of the DPD model 
is then like RD, which is solvable exactly, but with the 
additional ingredient of a density q of blocked sites. 

One can see that, within this approximation, the prob- 
ability of having a column with height h after N growth 
attempts is given by 
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where s = 1/L is the probability of attempting to 
growth a column and the usual approximation s r (l — 
s) N ~ r N\/[(N - r)lr\] « (Ns) r exp(-Ns)/r\ has been 
made. 

From the probability (0), one can calculate the inter- 
face width W 2 = (h 2 }~ (h) 2 and then the time derivative, 
which leading terms are 
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where t = Ns = N/L is the time in the units used in 
Ref. 0. This formula gives the exact time evolution of 
for any q (not only at q c = 0.539) and is valid for 
times (« 1. For times t > 1 differences between neigh- 
bouring columns are likely to be larger than 2 resulting 
in horizontal correlations and the break down of (|2j) . A 
comparison of Eq.(||) with numerical simulations of the 
DPD model is presented in Figure 1. 

Our calculation suggests that the exponential term in 
the ansatz of Ref. is actually produced by the usual 
random deposition-like dynamics, which occurs in any 
growth model Q for short times. 
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FIG. 1. Numerical results for the DPD model in a system 
of size L = 2 13 for q c — 0.539 (circles) and q = 0.3 (squares). 
Continuous lines correspond to Eq.(2) and fit the data for 
t < 1. For larger times our approximation is not valid any 
longer and the power law £ 2/3_1 takes over with (3 = 0.623 and 
P = 0.3 for q c — 0.539 and q = 0.3 respectively (dotted lines). 


